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Ferromagnetic order of nuclear spins coupled to conduction electrons: 
a combined effect of the electron-electron and spin-orbit interactions 
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We analyze the ordered state of nuclear spins embedded in an interacting two-dimensional electron 
gas (2DEG) with Rashba spin-orbit interaction (SOI). Stability of the ferromagnetic nuclear-spin 
phase is governed by nonanalytic dependences of the electron spin susceptibility x %3 on the mo- 
mentum (q) and on the SOI coupling constant (a). The uniform (q — 0) spin susceptibility is 
anisotropic (with the out-of-plane component, , being larger than the in-plane one, X° x , by a 
term proportional to J7 2 (2fci?)|a|, where U(q) is the electron-electron interaction). For q < 2m*\a\, 
corrections to the leading, U 2 (2kp)\a\, term scale linearly with q for x xx and are absent for \ zz ■ 
This anisotropy has important consequences for the ferromagnetic nuclear-spin phase: (i) the or- 
dered state-if achieved-is of an Ising type and (ii) the spin-wave dispersion is gapped at q = 0. 
To second order in (7(g), the dispersion a decreasing function of q, and anisotropy is not sufficient 
to stabilize long-range order. However, renormalization in the Cooper channel for q <C 2m*|a| is 
capable of reversing the sign of the (^-dependence of x xx and thus stabilizing the ordered state. We 
also show that a combination of the electron-electron and SO interactions leads to a new effect: 
long-wavelength Friedel oscillations in the spin (but not charge) electron density induced by local 
magnetic moments. The period of these oscillations is given by the SO length 7r/m*|a|. 

PACS numbers: 71.10.Ay,71.10. Pm, 75.40. Cx 
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I. INTRODUCTION 

Spontaneous nuclear spin polarization in semiconduc- 
tor heterostructures at finite but low temperatures has 
recently attracted a considerable attention both on the 
theoretical^— and experimental^ sides. Apart from 
a fundamental interest in the new type of a ferromag- 
netic phase transition, the interest is also motivated by 
an expectation that spontaneous polarization of nuclear 
spins should suppress decoherence in single-electron spin 
qubits caused by the hyperfme interaction with the sur- 
rounding nuclear spinsii 2 -, and ultimately facilitate quan- 
tum computing with single-electron spins £i£ 

Improvements in experimental techniques have lead 
to extending the longitudinal spin relaxation times in 
semiconductor quantum dots (QDs) to as long as ls£~— 
The decoherence time in single electron GaAs QDs has 
been reported to exceed 1/xs in experiments using spin- 
echo techniques at magnetic fields below lOOmT ; 11 ' 12 
whereas a dephasing time of GaAs electron-spin qubits 
coupled to a nuclear bath has lately been measured to 
be above 200^sJ^ Still, even state-of-the-art dynamical 
nuclear polarization methods^— allow for merely up to 
60% polarization of nuclear spins^ whereas polarization 
of above 99% is required in order to extend the electron 
spin decay time only by one order of magnitude Full 
magnetization of nuclear spins by virtue of a ferromag- 
netic nuclear spin phase transition (FNSPT), if achieved 
in practice, promises a drastic improvement over other 
decoherence reduction techniques. 

The main mechanism of the interaction between nu- 
clear spins in the presence of conduction electrons is the 



Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction^ 
The effective Hamiltonian of the RKKY interaction be- 
tween on-site nuclear spins of magnitude I 

ffRKKY - ~£ J^(T,T')r(r)P(T'), (1.1) 



is parameterized by an effective exchange coupling 

J^(r,r') = ^X lJ (r,r'), (1.2) 

where A is the hyperfme coupling constant, n s is the 
number density of nuclear spins, and 



X lJ (r, 



l/T 



dr(T T S l (r,r) S> (r',0)) (1.3) 



is the (static) correlation function of electron spins. 
[Hereafter, we will refer to x y (r, r')-and to its 
momentum-space Fourier transform-as to "spin suscep- 
tibility" , although it is to be understood that this quan- 
tity differs from the thermodynamic susceptibility, de- 
fined as a correlation function of electron magnetization, 
by a factor of /i 2 ^, where /is is the Bohr magneton.] It is 
worth emphasizing that x y (r, r') contains all the effects 
of the electron-electron interactio n 1 ' 2 this circumstance 
has two important consequences for the RKKY coupling. 
First, the electron-electron interaction increases the uni- 
form spin susceptibility which should lead to an enhance- 
ment of the critical temperature of the FNSPT, at least at 
the mean-field level. Second, stability of the nuclear-spin 
ferromagnetic order is controlled by the long- wavelength 



2 



behavior of the magnon dispersion w(q) which, in its 
turn, is determined by x lJ (q) at g 0. In a spin- 
isotropic and translationally invariant system, 



(1.4) 



with x lJ — $ijXi while the magnetization is given by 



M(T) = MJV J 



d D q 1 
, 16BZ (2tt)D e ^)/T _ I 



(1.5) 



where /ijv is the nuclear-spin magneton (we set fcs = h = 
1 throughout the paper). The second term in Eq. (fjjjj) 
describes a reduction in the magnetization due to ther- 
mally excited magnons. In a free two-dimensional elec- 
tron gas (2DEG), x(o) is constant for q < 2kp, and thus 
the magnon contribution to M(T) diverges in the 5^0 
limit, which means that long-range order (LRO) is un- 
stable. However, residual interactions among the Fermi- 
liquid quasiparticles lead to a non-analytic behavior of 
the spin-susceptibility: for q <C kp, x(<?) = x(0) + Cq, 
where both the magnitude and the sign of C depend on 
the strength of the electron-electron interaction ! 20 ' 21 In 
two opposite limits-at weak-coupling and near the Stoner 
instability^-the prefactor C is positive which, accord- 
ing to Eqs. (|1.4I) and (jl.5[) . means that LRO is unsta- 
ble. However, C is negative (and thus the integral in 
Eq. (|1.5|) is convergent) near a Kohn-Luttinger super- 
conducting instability^^ also, in a generic Fermi liquid 
with neither strong nor weak interactions C is likely to 
be negative due to higher-order scattering processes in 
the particle-hole channeli 2 ^— 

The spin-wave-theory argument presented above is 
supported by the analysis of the RKKY kernel in real 
space. A linear-in-g term in x(<?) corresponds to a dipole- 
dipole-like., 1/r 3 term in x( r ) (see Sec. Mil) . If C > 0, 
the dipole-dipole interaction is repulsive, and the ferro- 
magnetic ground state is unstable; vice versa, if C < 0, 
the dipole-dipole attraction stabilizes the ferromagnetic 
state. 

It is worth noting here that even finiteness of the 
magnon contribution to the magnetization does not guar- 
antee the existence of LRO. Although the Mermin- 
Wagner theorem 2 ^ in its original formulation is valid only 
for sufficiently short-range forces and thus not applicable 
to the RKKY interaction, it has recently been proven 2 ^ 
that magnetic LRO is impossible even for the RKKY in- 
teraction in D < 2. From the practical point of view, 
however, the absence of LRO in 2D is not really detri- 
mental for suppression of nuclear-spin induced decoher- 
ence. Indeed, nuclear spins need to be ordered within the 
size of the electron qubit (a double QD system formed 
by gating a 2DEG) as well as its immediate surround- 
ing such that there is no flow of magnetization. Since 
fluctuations grow only as a logarithm of the system size 
in 2D, it is always possible to achieve a quasi-LRO at 
low enough temperatures and on a scale smaller that the 



thermal correlation length. In addition, spin-orbit inter- 
action (SOI)-which is the main subject of this paper, see 
below-makes a long-range order possible even in 2D. 28 

The electron spin susceptibility in Eq. (11.41) was as- 
sumed to be at zero temperature. First, since the nu- 
clear spin temperature is finite, the system as a whole 
is not in equilibrium. However, a time scale associ- 
ated with 'equilibration' is sufficiently long to assume 
that there is no energy transfer from the nuclear- to 
electron-spin system. Second, if the electron tempera- 
ture is finite, the linear q scaling of x(<?) is cut off at 
the momentum of order T/vp = 1/Lt- For q <C 1/Lt, 
x{T,q) oc T + 0(v 2 F q 2 /T) such that u(q) oc q 2 and, 
according to Eq. (|1.5p . spin waves would destroy LRO. 
However, at low enough temperatures the thermal length 
Lt is much larger than a typical size of the electron qubit 
Lq. (For example, Lp ~ 1mm at T ~ lmK.) Therefore, 
q > 1/Lq ^> 1/Lt = T/vp and, indeed, the electron 
temperature can be assumed to be zero. 

In practically all nuclear-spin systems of current in- 
terest, such as GaAs or carbon-13 nanotubes, spin-orbit 
interaction (SOI) plays a vital role. The main focus of 
the paper is the combined effect of the electron-electron 
and SO interactions on the spin susceptibility of 2DEG 
and, in particular, on its q dependence, and thus on the 
existence/stability of the nuclear-spin ferromagnetic or- 
der. 

The interplay between the electron-electron and SOIs 
is of crucial importance here. Although the SOI breaks 
spin-rotational invariance and thus may be expected to 
result in an anisotropic spin response, this does not hap- 
pen for the Rashba and Dresselhaus SOIs alone: the 
spin susceptibility of free electrons is isotropic [up to 
cxp(— Ep/T) terms] as long as both spin-orbit-split sub- 
bands remain occupied^ The electron-electron interac- 
tion breaks isotropy, which can be proven within a Fermi- 
liquid formalism generalized for systems with SOI. 29 Spe- 
cific models adhere to this general statement. In partic- 
ular, x zz > X xx = X yy f° r a dense electron gas with the 
Coulomb interaction^ - 

In this paper, we analyze the q dependence of the spin 
susceptibility in the presence of the SOI. The natural 
momentum-space scale introduced by a (weak) Rashba 
SOI with coupling constant a (\a\ vp) is the difference 
of the Fermi momenta in two Rashba subbands: 



la 



2m* 



(1.6) 



where m* is the band mass of 2DEG. Accordingly, the de- 
pendence of x^ on <7 is different for q above and below q a ; 
in the latter case, it is also different for the out-of-plane 
and in-plane components. To second order in electron- 
electron interaction with potential U(q), the out-of-plane 
component is independent of q for q < q a : 

6 x zz (q,a) = 2xou 2 J-^. (1.7a) 
On the other hand, the in-plane component scales lin- 
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early with q even for q < q a 



S X xx (q,a) = 5 X yv (q,a) 



Xoul kf 



\a\kp 
3E F 



4 vpq 



(1.7b) 



In Eqs. (|1.7a|1.7bp . u q = m*U(q)/4:TT, kp is the Fermi 
momentum, Ep — k F /2m* is the Fermi energy, \o — 
m* /tt is the spin susceptibility of a free 2DEG, and Sx^ 
denotes a nonanalytic part of x IJ - For q a <C q <C hp, the 
spin susceptibility goes back to the result of Ref. [2CJ valid 
in the absence of the SOI: 



Sx lJ (q,a = 0) = %^Xou| fcF -|-. 



(1.8) 



Note that the subleading term in q in Eq. (|l.Tb[) differs 
by a factor of 2 /3 from the leading term in q in Eq. (|1.8|) . 
There is no contradiction, however, because Eqs. (jl.8[) 
and (|1.7b[) correspond to the regions of q < q a and q » 
q a , respectively. 

Equations (|1.7a[) and (jl.7bl) show that the uniform spin 
susceptibility is anisotropic: Sx zz (0,ce) = 2(5x^(0, a). 
This implies that the RKKY coupling is stronger if 
nuclear spins are aligned along the normal to the 
2DEG plane, and thus the nuclear-spin order is of the 
Ising type. In general, a 2D Heisenberg system with 
anisotropic exchange interaction is expected to have 
a finite-temperature phase transition^ In an anisotropic 
case, the dispersion of the out-of-plane spin-wave mode 
is given by^^ 



uj(q) 



if 

4n, 



(1.9) 



with q _L z. Ising- like anisotropy implies a finite gap in 
the magnon spectrum. In our case, however, the situa- 
tion is complicated by the positive slope of the linear q 
dependence of the second-order result for X xx (q)i which- 
according to Eq. (I1.9p -translates into to(q) decreasing 
with q. Combining the asymptotic forms of x lJ from 
Eqs. (|1.7all.7bj) . and (TO)) together, as shown in Fig. [Q 
we see that oj(q) is necessarily negative in the interval 
q a <C q -C kp, and thus LRO is unstable. Therefore, 
anisotropy alone is not sufficient to ensure the stabil- 
ity of LRO: in order to reverse the sign of the q de- 
pendence, one also needs to invoke other mechanisms, 
arising from higher orders in the electron-electron inter- 
action. We show that at least one of these mechanisms- 
renormalization in the Cooper channel-is still operational 
even for q <C q a and capable of reversing the sign of the 
g-dependence is the system is close to (but not neces- 
sarily in the immediate vicinity of) the Kohn-Luttingcr 
instability. 

We note that the dependences of Sx^ on q in the pres- 
ence of the SOI is similar to the dependences on the tem- 
perature and magnetic field^ presented below for coni- 
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FIG. 1: (color online): A normalized dispersion of the out- 
of-plane spin-wave mode uj(q) = cj(q)/L4 2 /xo/4n s ] as a func- 
tion of the momentum. To second order in interaction (lower 
curve) uo(q) is necessarily negative for m\a\ <C q <C fc_F, and 
thus LRO is unstable. Solid parts of the curves corresponds to 
actual calculations; dashed parts are interpolations between 
various asymptotic regimes. Renormalization effects in the 
Cooper channel reverse the slope of u(q) (upper curve) and 
stabilize LRO. 



pleteness: 

Sx zz (T,a) 
Sx zz {B z ,a) 

8 X xx {T,a) 
8x xx [B x ,a) 



2 Xou 2kf 



X0«2fc p 



\a\kp 

3E F 
' H k F 

3E F 
a\kp T 



0(T 3 ) 
0(Al) 



3Ep Ep 
\a\k F 16 |A 



O (T d 



3E F 3tt E F 



(1.10) 



Here, A l = gii B Bi/2 and T, A, < \a\kp. As 
Eqs. (|1.7all.7b|) and (|1.10[) demonstrate, while nonan- 
alytic scaling of 8x zz with all three variables (q, T, B) 
is cut off by the scale introduced by SOI, scaling of 5x zz 
continues below the SOI scale. This difference was shown 
in Ref. @ to arise from the differences in the dependence 
of the energies of particle-hole pairs with zero total mo- 
mentum on the magnetic field: while the energy of such 
a pair depends on the SO energy for B||z, this energy 
drops out for Biz. 

In addition to modifying the behavior of x %3 f° r <Z ^ 9a > 
SOI leads to a new type of the Kohn anomaly arising due 
to interband transitions: a nonanalyticity of x lJ (?> a ) 
at q = q a . The nonanalyticity is stronger in x zz than 
in x xx - 8x zz (q ~ q a ) oc (q - q a f /2 0(<7 - q a ) while 

8x xx {q ~ qa) « {q ~ qa)°^ 2 ©(<? — q a ), where Q(x) is the 
step-function. Consequently, the real-space RKKY in- 
teraction exhibits long- wavelength oscillations x zz ( r ) <x 
cos(q a r)/r 3 and (r) oc sin(q a r)/r 4 , in addition to con- 
ventional Friedel oscillations behaving as s'm(2kp r)/r 2 . 
It is worth noting that the long-wavelength Friedel oscil- 
lations occur only in the presence of both electron-electron 
and SO interactions. 

This paper is organized as follows. In Sec.|TT]we derive 
perturbatively the electron spin susceptibility of inter- 
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acting 2DEG with the SOI as a function of momentum; 
in particular, Sees. Ill Al III Dl outline the derivation of all 
relevant second-order diagrams, Sec. Ill El is devoted to 
Cooper renormalization of the second order result, and 
in Sec. Ill Fl we show that, in contrast to the spin sus- 
ceptibility, the charge susceptibility is analytic at small q 
(as it is also the case in the absence of SOI) . In Sec. IIIII 
we derive the real-space form the of the RKKY interac- 
tion and show that it exhibits long-wavelength oscilla- 
tions with period given by the SO length 2Tr/q a . Details 
of the calculations are delegated to Appendices 1X110] In 
particular, the free energy in the presence of the SOI 
is derived beyond the Random Phase Approximation in 
Appendix |DJ The summary and discussion of the main 
results are provided in Sec. IIVI 



"Cooper-channel" processes ; 20,24 in which two fermions 
with initial momenta k and p move in almost oppo- 
site directions, such that k ps —p. Backscattering pro- 
cesses are further subdivided into those with small mo- 
mentum transfer, such that (k, — k) — > (k, — k), and 
those with momentum transfers near 2k f, such that 
(k, — k) — > (— k, k). In the net result, all q = contri- 
butions cancel out and only 2kp contributions survive. 
We will show that this also the case in the presence of 
the SOI. In what follows, all "q — diagrams" are to be 
understood as the q = channel of the backscattering 
process. 



A. Diagram 1 



II. SPIN SUSCEPTIBILITY OF INTERACTING 
ELECTRON GAS 



1. General formulation 



Dynamics of a free electron in a two-dimensional elec- 
tron gas (2DEG) in the presence of the Rashba spin-orbit 
interaction (SOI) with a coupling strength a is described 
by the following Hamiltonian 



H = 



P 

2m* 



+ a(p x a y ~p y a x ), 



where p = (p x ,p v ) is the electron momentum of an elec- 
tron, and <r is a vector of Pauli matrices. The interaction 
between electrons will be treated perturbatively. For this 
purpose, we introduce a Green's function 



G(P) = 



1 



iujp — H — Ep 



= J2^s(p)9s(P) (2.2) 



with 



and 



n-(p) = g 



9s(P) 



1 + '-{p y <7 x -p x a y ) 
P 



sap 



(2.3) 



(2.4) 



where P = (oj p ,p) with u> p being a fermionic Matsubara 
frequency, e p = p 2 /2m* — Ep, and s = ±1 is a Rashba 
index. 

The nonanalytic part of a spin susceptibility tensor to 
second order in electron-electron interaction is given by 
seven linear response diagrams depicted in Figs.[2][7l Due 
to symmetry of the Rashba SOI, X* J (q) = X %i (o)$ij an d 

In the following subsections, we calculate all diagrams 
that contribute to non-analytic behavior of the out-of- 
plane, \ zz \ an d in-plane, \ xx — X yV i components of 
the spin susceptibility tensor for small external moment 
(q <C kp) and at T = 0. In the absence of SOI, the non- 
analytic contributions to the spin susceptibility from in- 
dividual diagrams are determined by "backscattering" or 



The first diagram is a self-energy insertion into the 
free-electron spin susceptibility, see Fig. [2j There are 
two contributions to the nonanalytic behavior: (£) from 
the region of small momentum transfers, i.e., q <C kp, 



(2.1) x^ = o(9)=2C/ 2 (0) 



Tr[G{P)G(P + Q)} 



Q JK JP 

x Tr[G(K + Q)a l G(K)G(K + Q)G(K)a j ] 

(2.5a) 

and (ii) from the region of momentum transfers close to 
2kp, i.e., |k — p| w 2kp and q <C kp, 

xl q =2k P (?) =2t/ 2 (2M / / [ Tv[G(K + Q)G(P + Q)] 
JqJkJp 

x Tt[G(K + Q)a l G{K)G{P)G{K)a J ]. 

(2.5b) 

Here, K = (wfc,k) and J„ = (2n)~ 3 J dujkd 2 k (and 
the same for other momenta). The time component of 
Q = (fl, q) is equal to zero throughout the paper. Since 
the calculation is performed at T = 0, there is no differ- 
ence between the fermionic and bosonic Matsubara fre- 
quencies. A factor of 2 appears because the self-energy 
can be inserted either into the upper or the lower arm 
of the free-electron susceptibility. As subsequent analy- 
sis will show, a typical value of the momentum transfer 
q is on the order of either the external momentum q or 
the "Rashba momentum" q a [cf. Eq. (11.61) ]. whichever is 
larger. In both cases, q <C kp while the momenta of 
both fermions are near kp, thus we neglect q in the an- 
gular dependencies of the Rashba vertices: fl s (k + q) ~ 
O s (k + q) « n s (k) = {l + s(sm9 k qcr x -cos6 k qa v )}/2 with 
9 a b = Z(a, b). [The origin of the x-axis is arbitrary and 
can be chosen along q.] Also, we impose the backscatter- 
ing correlation between the fermionic momenta: k = — p 
in the 2fcp-part of the diagram. With these simplifica- 



5 



tions, we obtain 

M=2^<0>/f/^/f^, 

x Iimnr(tt, 9kg,q, q)IL st (Cl, q), (2.6a) 

x Ii mnr (Q, 9 k g, q, g)n at (n, g), (2.6b) 
where summation over the Rashba indices is implied, 
a\l nnr ee Tr[Q^kyn m (k)Q n (k)O r (ky], (2.7a) 

b st ee Ir[fi s (p)O t (p)] = (1 + st)/2, (2.7b) 
a« sr ee Tr[^(k)a l r! m (k)^(-k)O r (k)^], (2.7c) 
b nt ee Tr[fi n (-p)O t (p)] = (1 - nt)/2 (2.7d) 




x .g;(w fc ,k + q)g m (u; fc ,k)g„(w fc + f2, k + q).g r (w fe , k) 

(2.7e) 

and, finally, the partial components of the particle-hole 
bubble are given by 

x g s (uj p , v)gt{u P + 0, p + q) 



27r y/vjrq 2 + (0 + i(i - s)ak F ) 2 ' 

For the derivation of the particle-hole bubble, see, e.g. 
Ref. S Calculation of other common integrals is pre- 
sented in Appendix lAl 

The main difference between the out-of-plane and in- 
plane components is in the structure of the "quaternion" , 
defined by Eq. (|2.7ep and calculated explicitly in Ap- 
pendix O [cf. Eq. (IA3|) ] . The dependence of Ii mnr on 
the external momentum q enters only in a combination 
with the SOI coupling as vf<] cos 9kq + (s — s')akp, where 
s, s' G {I, to, 7i, r}. Combinations of indices I, m, n, r are 
determined by the spin vertices a 1,1 and are, therefore, 
different for the out-of-plane and in-plane components. 
The out-of-plane component contains only such combina- 
tions {1, 77i, 7i, r} for which the coefficient s — s' is finite. 
Therefore, the SOI energy scale is always present and, 
for q <§C q a , one can expand in q/q a - The leading term in 
this expansion is proportional to |a| but any finite-order 
correction in q/q a vanishes. In fact, one can calculate 
the entire dependence of xf z on q (what is done in Ap- 
pendix [Bj and show that \f z is indeed independent of q 



P + QJ 




K + Q.l 

p + at 




K + Q, I 



FIG. 2: Diagram 1. Top: small-momentum transfer part. 
Bottom: 2fcF-momentum transfer part. K, s denotes a 
fermion from Rashba subband s — ±1 with "four-momentum" 
K=(u k ,k). 

for q < q a (and similar for the remaining diagrams). On 
the other hand, some quaternions, entering the in-plane 
component, have s = s' and thus do not contain the SOI, 
which means that one cannot expand in q/q a anymore. 
These quaternions provide linear-in-g dependence of X\ x 
even for q < q a , where the slope of this dependence is 
2/3 of that in the absence of the SOI. This is the origin 
of the difference in the q dependencies of x zz and \ xa \ as 
presented by Eqs. (|1.7al) and (ll.7b|) . 

The evaluation of the out-of-plane and in-plane part of 
diagram 1 is a subject of the next two subsections. 

2. Diagram 1: out-of-plane component 

We begin with the out-of-plane component of the spin 
susceptibility, in which case af^ nr = [1 + mr + n(m + 
r) — l{m + 7i + r + mnr)}/8 and af^j sr = [1 + 77ir — s(m + 
r) + l(s — to — ?' + T777\s)]/8. Summation over the Rashba 
indices yields 




x [(/+___ + /-+++)n 

+ /+-+-!!+_ + I_ + _ + n_+], (2.8b) 
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where IIo = n ++ = II 

As we explained in Sec. Ill A 11 the quaternions in 
Eqs. (|2.8a[) and (|2.8b[) contain q only in combination with 
q a . Therefore, for q <C q a , the leading term is obtained 
by simply setting q — 0, upon which the remaining inte- 
grals can be readily calculated. The results are given by 
Eqs. (TAT)) and (|A"8j) . so that 



■v zz 



2 \ a \ k F 

u Xo 



and 



Xl,q=2k f 



U 2k F X0 



3E F 

\a\kp 
3E F ' 



(2.9a) 



(2.9b) 




2a) 



In fact, it is shown in Appendix [B] that Eqs. (|2.9a[) and 
(|2.9b|) hold for any q < q a rather than only for q = 0. 



3. Diagram 1: in-plane component 

The in-plane component of the spin susceptibility dif- 
fers substantially from its out-of-plane counterpart due 

the angular dependence of the traces aj J r and a] 1 

which, for the in-plane case, read as 



2h) 



FIG. 3: Diagram 2. Top: small-momentum transfer part. 
Bottom: 2fc_F-momentum transfer part. 



Notice that the linear-in-g dependence survives in the 
in-plane component of the spin susceptibility even for 
<Z < Qa- Similar behavior was found in Ref. 13 for the 
temperature dependence of the uniform spin susceptibil- 
ity in the presence of the SOI. 



B. Diagram 2 



Iranr 



Diagram 2, shown in Fig. [5J is a vertex correction to 
le spin susceptibility. As in the previous case, there 
vo regions of momentum transfers relevant for the : 
(2.Mr)alytic behavior of the spin susceptibility: the q = 



= i [1 + mr + n{m + r) — l(m + n + r + mnr) cos 29 k ] 

1 the spin susceptibility. As in the previous case, there are 

= g [1 + nir - s(m + r) + l(s — m — r + mrs) cos 28 k ]^ w0 re gi ons f momentum transfers relevant for the non- 



(For the sake of convenience, we choose the x axis to 
be perpendicular to q when calculating all diagrams for 
X xx ■) Summing over the Rashba indices, one arrives at 



region, where 

X2U0 =u 2 (o) 



Xff«=o =4^(0) / ^ 



dfl f dOkq 



2?r 



qdq 
2tt 



Q JK 



Ti[G(P)G(P + Q)} 



[sin 2 e M (i+— +I-+++ )n 



x Tt[G(K + Q)G(K + Q + Q)<j l G{K + Q)G(K)a j ], 

(2.13a) 



and 



XX 

Xl,q=2k F 



COS 2 6 k q I— 
dOkq 



--2U 2 (2k F ) 



dfl 
2^ 



2vr 



.)n ] (2.11a) 

qdq 
2vr 



and the 2fcp-region, where 



x [sin 2 9 k q(I+ +/- +++ )n 



COS 2 6kq(I-i 



)n 



X2,q=2k f 



Explicitly, 



U 2 (2k F ) [ [ [ Tt[G(K + Q)G(P + Q)} 

JQ JK JP 



x Tr[G(K + Q)G(P + Q)a l G(P)G(K)a ] }. 

(2.13b) 



+ sin 2 6» fe g(7 + _ + _n+_ + 

+ cos 2 (9 &5 (i ++ _ + n_ + + i__ + _n+_)]. 

(2.11b) 

Details of the calculation are given in Appendix. IA 2 bl 
here we present only the results in the interval q < q a : 



X2U =^(0) / ^ 



dn r d0 k , 



kq 



2n 



qdq 
2tt ' 



•V I 

Iranr 



x J l in n r (n,9 k q, q ,q)nst(n, q ), (2.i4a) 



XX 

Xl,q=0 



XX 

Xl,q=2k f 



4-„ 2 v 2 VF<i 
2 Xl 'i= + UoXo 9^W 



Xiq = 2k F = U\2k F ) I - 



dn r dd 



'kq 



2tt 



2tt 



\a\ki 



6E F 



2 v F q 
9tt 



2 Al,<?= 2kp 

ulk F Xo 



2 2 v F q 

\a\k F 2 v F q\ 
QE F 9?r E f ) ' 



(2.12a) 



(2.12b) 



where 



X3 
Irsm 



X hrnni^l-, dkq-, q, q)I rS t{^, Okq, <7, -<?), 

(2.14b) 



= Tr[0 ; (k)an(k)CT 4 0„(k)^ r (k)CT J ], (2.15a) 
Tr[Oi(k)O r .(-k)o- i O s (-k)O m (k)(r i ] ) (2.15b) 
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de. 



kq 



2tt 



Jlmnr(Cl, 9 k q, <?, O) — , , , 

Z7T J Zn 

x 5r ; (w fc + Cl, k + q)g m (wfc + Cl, k + q + q) 

x i g„(w fe + fi,k + q) i g r (wfc,k), (2.15c) 



d0. 



kq 



du>„ f de 



ZTT J ZTT J ZTT 

x gi(uj k ,k + q)g m (uj k ,k)g n (uj k + ft,k + q), (2.15d) 

As before, summation over the Rashba is implied. Inte- 
grals (|2.15c[) and (|2. 15d[) are derived in Appendix [A] 

Traces entering the q = part of the out-of-plane and 
in-plane components are evaluated as 

zz 1 + nr — m(n + r) + l(m — n — r + mnr) 



XX 



(1 + lm)(l + nr) + (I + m)(n + r) cos29 k q 



(2.16) 



Summing over the Rashba indices and using the symme- 
try properties of hmnr and Jimnr, it can be shown that 
the q — parts of diagrams 1 and 2 cancel each other 



ij ij 

X2,q=0 — ~ Xl,g=0l 



(2.17) 



which is also the case in the absence of the SOI. 20 There- 
fore, we only need to calculate the 2fcp-part of diagram 2. 



1. Diagram 2: out-of-plane component 

Summation over the Rashba indices with the coefficient 
cf* sm = [1 + mr — s(m + r) + l(s — r — m + mrs)]/8 for 
the out-of-plane part gives 

_ TT 2 (0h v f ^ f d0 kg f qdq 

X2,q=2k F - U (2k F )j — J — J — 

x Okq, q, q)I-+- (Cl, 9kq, q, -q) 

+ j+__(n, e k q, q, q)i-++{n, e kq , q , -§) + (q -> -§)], 

(2.18) 

where (q — > — g) stands for the preceding terms with 
an opposite sign of momentum. Integrating over q and 
Cl at q = 0, yields [cf. Eq. (|A9j) ]. 



X2.q=2k F — u 2k F XCf 



\a\kp 
3E F ' 



(2.19) 



Again, an exact calculation at finite g proves that this 
results holds for any q < q a . 



2. Diagram 2: in-plane component 

The in-plane component comes with a Rashba coeffi- 
cient cf^ sr = [(1 - lr)(l - ms)(l - r)(m - s) cos 29 kq ]/8, 



K, rn 




3) 



1) 



FIG. 4: Top: diagram 3. Bottom: diagram 4. The momen- 
tum transfer q in both diagrams can be either small or close 
to 2k F . 



such that 



XX 

A2,<j=2fc 



U 2 (2k F ) 



dci r de 



'kq 



qdq 
2vr 



2tt J 2tt 

x {sin 2 6 kq \I + ^ + (Cl, 6 k q, q, q)I_ + _((l, fe5 , q, -q) 
+ 1+ — (O) ^ > <7, <D-T-++(^> <z, -g)] 

+ cos 2 9 kq [I +++ (Cl, 9 kq , q, q)I (Cl, 9 kq , q, -q) 

+ (Cl, 9 M , q, q)I—+(Cl, d k q, q, -q)] 

+ ($->-£)}. (2.20) 

The first part, proportional to sin 2 9 kq -, contains the SOI 
coupling a. In this part, q can be set to zero, and the 
resulting linear-in-|a| part equals half of that for the out- 
of-plane component due to the integral over sin 2 9 kq -. On 
the other hand, in the term proportional to cos 2 9 kq -, the 
dependence on \a\ drops out upon integration over q, and 
the final result for q < q a reads as [cf. see Eq. (|A12|) ] 



XX 

X2,q=2kf 



',X2*q=2k F + u 2k F X0 



u 2k F X0 



2 v F q 

2 v F q\ 
6E F ' 9tt E F ) ' 



\OL\Kf 



+ 



(2.21) 



C. Diagrams 3 and 4 

We now turn to " Aslamazov-Larkin" diagrams, Fig. SI 
which represent interaction via fluctuational particle-hole 
pairs. Without SOI, these diagrams are identically equal 
to zero because the spin vertices are averaged indepen- 
dently and thus vanish. With SOI, this argument does 
not hold because the Green's functions now also contain 
Pauli matrices and, in general, diagrams 3 and 4 do not 
vanish identically. Nevertheless, we show here the non- 
analytic parts of diagrams 2 and 3 are still equal to zero. 

Diagrams 3 and 4 correspond to the following analyti- 
cal expressions: 

*3 J =/ / I U 2 (H)Tr[G{P-Q)G(P-Q)G(P)o-*] 
Jq J k Jp 

xTt[G(K + Q)G(K + Q)G(K)a j ], 

(2.22a) 
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X 3 = 



III U 2 {\q\)Tr[G{P)G{P + Q)G(P + Q)a l 

JQ JK J P 



xTr[G{K + Q)G(K + Q)G{K)<j j }. 

(2.22b) 

Note that the second trace is the same in both diagrams. 
In what follows, we prove that 

\V " (2-23) 
for both small and large momentum transfer q. 



The prime over / denotes that integration over the angle 
6kq is not yet performed as compared to Ii mn (fl, 6kq, q, q) 
defined by Eq. (plSd]) . 

Summing over the Rashba indices, one finds 



^df nm li mn (n,e kq ,q,q)=0 (2.28) 



1. Diagrams 3 and 4: out-of-plane components 



The out-of-plane case is straightforward. Evaluating 
the second traces in Eqs. (|2.22a|) and (|2.22b[) . one finds 
that they vanish: 

d\ nm = Tr[^(k)tt n (k)r! m (k)a z ] = 0, (2.24) 

for the q = case, and 

df nm = Tr[fi i (k)O n (-k)0 TO (k)cr 2 ] = 0, (2.25) 

for the q = 2kp case. Therefore, xf 2 = x\ z = 0- 

2. Diagrams 3 and 4: in-plane components 



and, therefore, the in-plane component at small momen- 
tum transfer vanishes. 

^(Iqllj^^^^ - pi: 



K, m 




5) 



FIG. 5: Diagram 5. The momentum transfer q is close to zero 
and |k — p| = 2ftF- 



For the in-plane part of the spin susceptibility, the 
proof is more complicated as the traces do not vanish 
on their own. To calculate the q = part, we need the 
following two objects 



df nm ee Tr[^(k)tt„(k)tt m (kK] 

= cos9kq(l + m + n + lmn)/A 



(2.26) 



and 



l'lmn{^ e kq,q,q) = 1^ J d ^k J de k gi(uj k , k + q) 
x g m (u k , k)g n (ujk + 0,k + q) 

zf2 — VFqcosOkq + VFqcosOkq + (I — n)akp 

x - . (2.27) 

iQ — VFqcosOkq + (m — n)akp 



The trace for the q = 2kp case turns out to be the 
same as for the q = one 



df nm = Tr[n i (k)fi„(-k)O m (k) C 7 a: ] = df n 



(2.29) 



However, in order to see the vanishing of the 2k f part, 
the integral over has to be evaluated explicitly with 
q = 2kp, i.e., 



Tn* f 

^Lni^^kq.q = 2k F ,q) = — / de k gi (uk , k + q)g m (u k , k)g n {uj k + 0, k + q) 

im*[i-e(u k )-e(u k + n)] 

[i(2wfe + Q) — vpq cos Okq — vpq — 2vpkp cos 9k q — (m + n)akp] [i(2uik + f2) — vpq — 2vpkp cos 9k q — {I + n)akp] ' 

(2.30) 
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^(ik- P i: 



a s 6) 



K + Q, I 



P-Q,t 



FIG. 6: Diagram 6. The momentum transfer q is close to zero 
and |k — p| = 2k F - 



where we used an expansion of £k+q around q = 2k F : 
£k+q ~ — tk + vp{q — 2k F ) + 2v F k F cos9 kq . Summing 
over the Rashba indices, we obtain 




a' 7a) 



76) 



K + Q,l 



P + Q, 



(2.31) 



Iran 



and, therefore, the 2k F part of the in-plane components 
of diagrams 3 and 4 is also equal to zero. 



D. Remaining diagrams and the final result for the 
spin susceptibility 

The remaining diagrams can be expressed in terms of 
the diagrams we have already calculated. 
Diagram 5 in Fig. [5] reads as 

X% =-AU{G)U{2k F ) [ [ f Tv[G{K + Q)a l G{K) 

JQ JK JP 

x G(K + Q)C 
= -4U(0)U(2k F 



FIG. 7: Diagram 5a (upper figure) and diagram 5b (lower 
figure). The transferred momenta are q = and k — p| = 2k F - 



are related as well. Explicitly, diagram 6 reads as 
Xl j = - 2U(0)U(2k F ) f f ( Tr[G(K + Q)a l G(K) 

JQ JK JP 

x G(K + Q)G{P - Q)o 3 G{P)G(P - Q)} 



with 



2tt J 2tt J 2tt 
x Ii mn (Q, 9 k q, q, q)I rst (-Q, 9 kq , -q, q) (2.35) 



x G(K + Q)G(P + Q)G{P)G{K)a ] 
dfl f dOkq 



2ir / 2tt 



qdq fi j 

n„- <* Imntsr 



Ilmnr^st 

(2.32) 



with 



ftLtsr = Tr[n ; (k)^fi m (k)Q n (k) 



(2.33) 



x n r (-k)o- J 'O s (-k)O t (-k)]. (2.36) 
On the other hand, for diagram 7b we obtain 

X% = - 2U(0)U(2k F ) III Tr[G{K + Q)a*G(K) 

JQ JK JP 

x G(K + Q)G(P + Q)G{P)a j G{P + Q)] 
dQ. f d9 kq f gdq-jj 

2yj- Imntsr 

(2.37) 



= - 2U(0)U(2k F ) / — 



2tt J 2ir 

X /j m n(fJ, 8kq, q, q)Irst (H, 9kq, S) 



with 



and q \k — p\ = 2k F . A factor of 4 appears because 
the "sunrise" self-energy can be inserted into either the 
lower or the upper arm of the bubble while each of the 
interaction line can carry momentum of either q = 
or q = 2k F . A minus sign is due to an odd number 
of fermionic loops. Upon summation over the Rashba 
indices, we obtain 



Imntsr 



Tr[O i (k)<r l fi m (k)O n (k) 

x n t (-k)Q s (-k)a j fl r (-k)}. (2.38) 
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x 5 



Xl,q=0 



In both cases, q <C |k — p| = 2k F . Using the symmetry 
propertyI rst (-fi, fc §, -q, -q) = -I_ r _ s _ 4 (f2, 0^, g, g) in 
xi j summing over the Rashba indices, and noticing that 
I +++ (ft,9 k q,q,q) = I (tt,9 k q,q,q), we arrive at 



(2.34) 



U(0)U(2k F ) U 2 {0) ' 
Diagrams 6 and 7b in Figs. [5] and [7J correspondingly, 



%i ii 
X 6 = Xn- 



(2.39) 



Finally, diagram 7a shown in Fig. [7J is related to dia- 
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gram 2 at small momentum transfer. Indeed, 



11 

X7a 



-2U(0)U(2k F ) / / / Tr[G(A + Q + Q)a l G(K 

JQ JK Jp 

x G{P)G(P + Q)G(K)a j G(K + Q)] 
d£L f d9kq f qdq 
2t7 _ 



contains only one energy scale. In other words, the \a\ 
term is renormalized via In \ a\ while the q is renormal- 
(reed via In q. For example, the third-order result for the 
2kf part of diagram 1 (Fig. [2J reads as (for simplicity, 
we assume here a contact interaction with U (q) = const) 



2U(0)U(2k F ) 



2n 



2tt 



-h tj 



(2.40) 



with 



hnnstr 



Tr[O i (k)Q ro (k)<r i O n (k) 



(2.41) 



where again q <C |k — p| = 2k p. After summation over 
the Rashba indices, this diagram proves related to the 
small-momentum part of diagram 2 as 



X5a 



X2,q=0 



U(0)U(2k F ) U 2 (0) 



(2.42) 



The results of this section along with Eq. (|2.17l) show 
that the sum of all diagrams proportional to U(0)U(2kp) 
cancel each other 



X 5 + Xi + X7a + X 7b = 



(2.43) 



Therefore, as in the absence of SOI, the non-analytic part 
of the spin susceptibility is determined only by the Kohn 
anomaly at q = 2k p. 

Summing up the contributions from diagrams 1 — 3, we 
obtain the results presented in Eqs. (|1.7ap and (|1.7b[) . 



E. Cooper-channel renormalization to higher 
orders in the electron-electron interaction 

An important question is how the second-order results, 
obtained earlier in this Section, are modified by higher- 
order effects. In the absence of SOI, the most impor- 
tant effect-at least within the weak-coupling approach- 
is logarithmic renormalization of the second-order re- 
sult by to the interaction in the Cooper channel. As 
it was shown in Ref. 0, this effect reverses the sign of the 
q dependence due to proximity to the Kohn-Luttinger 
superconducting instability; the sign reversal occurs at 
q = e Trl/vf R3 7ATkl/vf, where Tkl is the Kohn- 
Luttinger critical temperature. For momenta below the 
SO scale (q a ), X zz ceases to depend on q but \ xx still 
scales linearly with q. What is necessary to understand 
now is whether the linear-in-g term in \ xx renormalized 
in the Cooper channel. The answer to this question is 
quite natural. The \a\- and q terms in the second-order 
result for X xx [Eq. ((L7b|) ] come from different parts of 
diagram: the \a\ term comes from q independent part 
and vice versa. Starting from the third order and be- 
yond, these two terms acquire logarithmic renormaliza- 
tions but the main logarithm of these renormalizations 



xZ 



i2xo 



=2k f 



\a\kf 



In- 



A 



2 v F q j A 
37r Ep VFq 



Ep \a\kp 

(2.44) 

where u = m*U/Air and A is the ultraviolet cutoff. De- 
tails of this calculation are given in Appendix O It is 
clear already from this result the logarithmic renormal- 
ization of the q term in x xx remains operational even for 
q < q a , with consequences similar to those in Ref. [H 



F. Charge susceptibility 

In the absence of SOI, non-analytic behavior as a func- 
tion of external parameters-q, T, H-is present only in 
the spin but not charge susceptibility 20 i 33 i 34 An inter- 
esting question is whether the charge susceptibility also 
becomes non-analytic in the presence of SOI. We answer 
this question in the negative: the charge susceptibility 
remains analytic. To show this, we consider all seven 
diagrams replacing both spin vertices by unities. The 
calculation goes along the same lines as before, thereby 
we only list the results for specific diagrams; for q <C q a , 



Sxl = 


-&xl = 


xo , 2 , 2 \ v Fq 




-5x1 = 


xo i 2 2\ v pq 


5x1 = 


-5x1 = 


xo VFq 

-^ U ° U2k *Ep-' 



(2.45) 



whereas \7 = on its own {Xra = —Xn)- First, we im- 
mediately notice that SOI drops out from every diagram 
even in the limit q <C q a - Second, the sum of the non- 
analytic parts of all the charge susceptibility diagrams is 
zero, S\ c = 0, as in the case of no SOI. 



III. 



RKKY INTERACTION IN REAL SPACE 



A nonanalytic behavior of the spin susceptibility in 
the momentum space leads to a power-law decrease of 
the RKKY interaction with distance. In this Section, 
we discuss the relation between various nonanalyticities 
in x lJ ( ( ?) an d the real-space behavior of the RKKY in- 
teraction. We show that, in addition to conventional 
2kF Friedel oscillations, a combination of the electron- 
electron and SO interactions lead to a new effect: long- 
range Friedel-like oscillations with the period given by 
the SO length. 
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A. No spin-orbit interaction 



First, we discuss the case of no SOI, when the spin 
susceptibility is isotropic: X y (?) = $ijX(o)- T° r f ree 
electrons, the only non-analyticity in Xo(q) is the Kohn 
anomaly at q = 2k p, which translates into Friedel oscil- 
lations of the RKKY kernel; in 2D, and for k F r > l, 35 



Xo(r) 



Xq sin (2k F r) 
2ir r 2 



(3.1) 



One effect of the electron-electron interaction is a loga- 
rithmic amplification of the Kohn anomaly (which also 
becomes symmetric about the q = 2kp point): x{q ~ 
2kp) oc y/\q — 2kp\ In \q — 2kp\£& Consequently, x( r ) 1S 
also enhanced by logarithmic factor compared to the free- 
electron case: x( r ) ^ sm(2kpr)\n(kpr)/r 2 . 

Another effect is related to the nonanalyticity at small 
T- x(o) = Xo + CqM To second order in the electron- 
electron interaction [cf. Eq. 



C 2 



4xo 2 
-u 



3irki 



2k F i 



(3.2) 



however, as we explained in Sec. HI both the magnitude 



and sign of C can changed due to higher-order effects. 
(Cooper channel renormalization leads also to multiplica- 
tive In q corrections to the linear-in-g term; those corre- 
spond to multiplicative lnr renormalization of the real- 
space result and are ignored here.) 

In 2D, x( r ) is related to x(<?) via 



x(r) 



1 

2^ 



dqqx(q)Jo(qr). 



(3.3) 



Power-counting suggests that the q term in x(o) trans- 
lates into a dipole-dipole-like, 1/r 3 term in x( r )- To see 
if this indeed the case, we calculate the integral 



A 



dqq Jo(qr) 



(3.4) 



with an arbitrary cutoff A, and search for a universal, 
A-independent term in the result. If such a term exists, 
it corresponds to a long-range component of the RKKY 
interaction. Using an identity xJq(x) = j^{xJi{x)) and 
integrating by parts, we obtain 



(Ar) 2 Ji(Ar) 



Ar 



dxxJi(x) 



(Ar)Vi(Ar) 



7rAr 



r 



{Ji(Ar)H (Ar) - Jo(Ar)Hi(Ar)} 



(3.5) 



where H„(n) is the Struve function. The asymptotic ex- 
pansion of the last term in the preceding equation indeed 
contains a universal term 



7rAr 



{J!(Ar)H (Ar) - ^(A^H^Ar)} 



l 



(3.6) 

where . . . stands for non-universal terms. A correspond- 
ing term in x( r ) reads 



xW 



c 

2nr 3 



(3.7) 



As a check, we also calculate the Fourier transform 
of the g-independent term in x lJ ■ The corresponding 
integral 



.4 



A 

dqqJo(qr) = — Ji(Ar 



(3.8) 



does not contain a A-independent term and, therefore, 
a constant term in x(q) does not produce a long-range 
component of the RKKY interaction, which is indeed the 
case for free electrons. 

Equation (|3.7[) describes a dipole-dipole-like part of 
the RKKY interaction that falls off faster than Friedel 



oscillations but is not oscillatory. [Incidentally, it is the 
same behavior as that of a screened Coulomb potential 
in 2D, which also has a q nonanalyticity at small q£^] 

In a translationally invariant system, -ffRKKY = 
-£t Err' X(r - r')lili,. Therefore, if C> 0, i.e., X (q) 
increases with q, the dipole-dipole interaction is repulsive 
for parallel nuclear spins and attractive for antiparallel 
ones. Since the 1/r 3 behavior sets in only at large dis- 
tances, the resulting phase is a helimagnet rather than an 
antiferromagnet. Vice versa, if C < 0, the dipole-dipole 
interaction is attractive for parallel spins. This corre- 
sponds precisely to the conclusions drawn from the spin- 
wave theory: a stable FM phase requires that v (q) > 0, 
which is the case if C < 0. 



B. With spin-orbit interaction 

C. Free electrons 

In a free electron system, the SOI splits the Fermi sur- 
face into two surfaces corresponding to two branches of 
the Rashba spectrum with opposite helicitics. Conse- 
quently, both components of the spin susceptibility in 
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the momentum space have two Kohn anomalies located 
at momenta 2k p = 2k p =p q a with q a — 2m* \a\ . To see 
this explicitly, we evaluate the diagonal components of 
X 13 (q) for q « 2k F 

Xo(q) = -Y] I |(k,s| ( r' i |k + q,i)| 2 .g t (w,k + q) ffs (w,k). 

s,t J K 

(3.9) 

For q « 2/cp, the matrix elements of the spin operators 
in the helical basis reduce to 



(k + q,t|^|k,s) 



st). 



Therefore, x" (?) contains only contributions from intra- 
band transitions 



X X o X (q) 



x zz (q) = - g+(w,k + q). g+ (w,k) 

JK 

ff_(w,k + q) S _(w,k). (3.11) 



K 



Each of the two terms in Eq. (13.111) has its own Kohn 
anomaly at q = 2k s Fl s = ±. In real space, this cor- 
responds to beating of Friedel oscillations with a period 
2n/q a . 

This behavior needs to be contrasted with that of 
Friedel oscillations in the charge susceptibility, where-to 
leading order in a-the Kohn anomaly is present only at 
twice the Fermi momenta in the absence of SOIi2& Con- 
sequently, the period of Friedel oscillations is the same 
as in the absence of SOI. (Beating occurs in the presence 
of both Rashba and Dresselhaus interactions^) This is 
so because, for q near 2kp , the matrix element entering 
X c (q) reduces to 



Kk + q,t|k,s) 



1 



(1 - st) . 



which implies that x° contains only contributions from 
intcrband transitions: 



Xl (q) 



2 f ff+ (w,k + q)0_(w,k). (3.12) 

JK 



The Kohn anomaly in Xo corresponds to the nesting con- 
dition e^+q = — e k ' which is satisfied only for q = 2k F . 



1. Interacting electrons 

The electron-electron interaction is expected to affect 
the 2/c^-Kohn anomalies in x xx and X zz m a wa Y similar 
to that in the absence of SOI. However, a combination of 
the electron-electron and SO interaction leads to a new 
effect: a Kohn anomaly at the momentum q a <C 2kp. 
Consequently, the RKKY interaction contains a compo- 
nent which oscillates with a long period given by the SO 
length Xso = 27r/g Q rather than the half of the Fermi 
wavelength. 



To second order in the electron-electron interaction, 
the full dependence of the electron spin susceptibility on 
the momentum is shown in App. |B]to be given by 



Sx xx (q)^ + ^Ke 



q a . q a 
— arcsm — 

q q 



1- 



Sx zz (q) =C 2 qRe 



1- 



q a . q a 

— arcsm — 



(3.13a) 



(3.13b) 



Equations (|3.13a[) and (|3.13b[) are valid for an arbitrary 
value of the ratio q/q a (but for q <C kp). For q 3> q a , 
both 8x xx and 6x zz scale as q. For q -C q a , o~x xx contin- 
ues to scale as q (but with a smaller slope compared to 
the opposite case), while 5x zz is q independent. The 
crossover between the two regimes is not continuous, 
however: certain derivatives of both 8x xx and 8x zz di- 
verge at q = q a . Expanding around the singularity at 
q = q a , one finds 



Sx xx =-^f + f 



Q(q a - q) 



e(«- fe )(i + f(l-i 

V 5 \q a 



s x z 



--c 2 



5/2 



0(q a - q) + &(q - q a ) ( 1 + b[ - 1 

Qa 



(3.14a) 

3/2 



(3.14b) 



where Q(x) is the step-function, C 2 = nC 2 q a /2 and b = 
4\/2/37r. The q dependencies of 5x xx and Sx zz are shown 
in Fig. |U 

The singularity is stronger in Sx zz <x (<7 — <Za) 3 ^ 2 whose 
second derivative diverges at q = q a , whereas it is only 
third derivative of Sx xx oc (q— q a ) 5 ^ 2 that diverges at this 
point. Both divergences are weaker than the free-electron 
Kohn anomaly \ oc (q — 2kp) 1 / 2 . 

We now derive the real-space form of the RRKY in- 
teraction, starting from X zz (r). Substituting Eq. (|3.14bj) 
into Eq. Q3.3P and noting that only the part proportional 
to (q/q a — l) 3 ^ 2 contributes, we arrive at the following 
integral 



X" 



Z7T 



dqqJo(qr) 1 

\q a 



3/2 



(3.15) 



where A is an arbitrarily chosen cutoff which does affect 
the long-range behavior of x zz ( r )- Replacing Jo (a;) by its 
large-x asymptotic form and q by q a in all non-singular 
and non-oscillatory parts of the integrand, we simplify 
the previous expression to 




(q 



q a )r-^ 



(3.16) 
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FIG. 8: (color online) The nonanalytic part of the electron 
spin susceptibility in units of (2/3ir)u2 kF (\ol\/vf)xo as a func- 
tion of the momentum in units of q a = 2m*|«|. i = x,z. 
Solid: exact results (|3.13ap and (|3.13b|) . Dashed: approxi- 
mate results (|3.14a[) and (|3.14b|) valid near the singularity at 
q = q a - 



Finally, the conventional, 2kp Friedel oscillations should 
be added to Eqs. Q3.17P and Q3.19P to get a complete r 
dependence. The dipole-dipole part and long-wavelength 
Friedel oscillations fall off faster then conventional Friedel 
oscillations. In order to extract the long-wavelength part 
from the data, one needs to average the measured \ l ° (r) 
over many Fermi wavelengths. Recently, 2k p oscillations 
in the RKKY interaction between magnetic adatoms on 
metallic surfaces have been observed directly via scan- 
ning tunneling microscopy. 40 Hopefully, improvements in 
spatial resolution would allow for an experimental verifi- 
cation of our prediction for the long-wavelength compo- 
nent of the RKKY interaction. 

As a final remark, we showed in Sec. Ill Fl that the 
charge susceptibility does not exhibit small-q nonanalyt- 
icities. This result also implies that the long- wavelength 
oscillations are absent in the charge susceptibility; there- 
fore, Friedel oscillations produced by non-magnetic im- 
purities contains only a conventional, 2kp component. 



Integrating by parts twice and dropping the high-energy 
contribution, we arrive at an integral that converges at 
the upper limit. The final results reads 



X ZZ (r) 



■Xo 



2 u\ kp cos (q a r) 



37T 2 kf 



(3.17) 



Equation (|3.17p describes long-wavelength Friedel-like 
oscillations which fall off with r faster than the usual 2kp 
oscillations. Notice that Eq. (|3.17p . while valid formally 
only for q a r 3> 1, reproduces correctly the dipole-dipole 
term [Eq. (|3.7p with C — C 2 ] in the opposite limit of 
q a r <C 1. Therefore, Eq. (|3.17l) can be used an extrap- 
olation formula applicable for any value of q a r. 

In addition to the Kohn anomaly at q = q a , the in- 
plane component also contains a non-oscillatory but non- 
analytic term, proportional to q. As it was also the case 
in the absence of SOI, this term translates into a dipole- 
dipole part of the RKKY interaction. Analysis of Sec. Mil 
fully applies here: we just need to replace the prefactor C 
in Eq. ([37FJ by 2C 2 /3, where C 2 is defined by Eq. (j3~2i 
The role of the cutoff A in Eq. (|3.4I) is now being played 
by q a , therefore, C — > 2C 2 /3 for r > q~ l . For r <C g" 1 , 
the prefactor is the same as in the absence of SOI. Sum- 
marizing, the dipole-dipole part of the in-plane RKKY 
interaction is 



2 2 f 1/r 3 , for q a r < 1 



3jT 2 u 2k F Xo x I 2/ir ^ imqaj , > ^3.18) 



The oscillatory part of x xx { r ) is obtained by the same 
method as for x zz ( r )i one only needs to integrate by 
parts three times in order to obtain a convergent inte- 
gral. Consequently, x xx ( r ) f aus off with r as 1/r 4 . The 
r-dependence of x xx ( r )> resulting from the SOI, is given 
by a sum of the non-oscillatory and oscillatory parts 



X XX (r)=xT-d(r)+ X o 



1 u\ kp sm(q a r) 
3tt 2 q a kp r 4 



(3.19) 



IV. SUMMARY AND DISCUSSION 

We have studied the nonanalytic behavior of the 
electron spin susceptibility of a two-dimensional elec- 
tron gas (2DEG) with SOI as a function of momentum 
q = |q| in the context of a ferromagnetic nuclear-spin 
phase transition (FNSPT). Similarly to the dependence 
on temperature and magnetic-field^ 4 - the combined ef- 
fect of the electro-electron and spin-orbit interactions 
affects two distinct components of the spin susceptibil- 
ity tensor differently. For q < 2m*|a|, where m* is 
the effective electron mass and a is the spin-orbit cou- 
pling, the out-of-plane component of the spin suscep- 
tibility, x zz (g, a), does not depend on momentum (in 
other words, momentum-dependence is cut off by the 
SOI), [cf. Eq. (|1.7al) ]. whereas its in-plane counterparts, 
X xx (q,cn) = x vv (<ii a )i scale linearly with q even below 
the energy scale given by the SOI [cf. Eq. (|1.7b[> ]. No- 
tably, both results are exact for q < 2m*\a\. 

Beyond second order in electron-electron interaction 
renormalization effects in the Cooper channel, being the 
most relevant channel in the weak coupling regime, start 
to play a dominant role. As we have shown in Sec. Ill El 
the leading linear-in-|a| term becomes renormalized by 
In I a I, while the subleading linear-in-g term acquires ad- 
ditional lng dependence. This behavior is a natural con- 
sequence of the separation of energy scales in each of 
the diagrams and suggests that, in general, X °c 
U n E l Ei ln"~ 2 E h where E { stands for a generic energy 
scale (in our case Ei = {\a\kp,vpq} but temperature or 
the magnetic field could be included as well). 

Our analysis of the spin susceptibility gives impor- 
tant insights into the nature of a FNSPT. First, the 
SOI-induccd anisotropy of the spin susceptibility implies 
that the ordered phase is of an Ising type with nuclear 
spins aligned along the z-axis since x zz > X xx ■ Second, 
the ferromagnetic phase cannot be stable as long as the 
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higher-order effects of the electron-electron interaction 
are not taken into account. In this paper, we focused only 
on one type of those effects, i.e., renormalization in the 
Cooper channel. Without Cooper renormalization, the 
slope of the magnon dispersion is negative, even though 
the magnon spectrum is gapped at zero-momentum, cf. 
Fig. [TJ This implies that spin-wave excitations destroy 
the ferromagnetic order. Only inclusion of higher-order 
processes in the Cooper channel, similarly to the mech- 
anism proposed in Ref. 1, leads to the reversal of the 
slope of the spin susceptibility in the (not necessarily im- 
mediate) vicinity of the Kohn-Luttinger instability, and 
allows for the spin- wave dispersion to become positive at 
all values of the momentum. This ensures stability of the 
ordered phase at sufficiently low temperatures^^ 

We have also shown that a combination of the electron- 
electron and SO interactions leads to a new effect: a 
Kohn anomaly at the momentum splitting of the two 
Rashba subbands. Consequently, the real-space RKKY 
interaction has a long- wavelength component with a pe- 
riod determined by the SO rather than the Fermi wave- 
length. 

Another issue is whether the SOI modifies the behav- 
ior of the charge susceptibility which is known to be ana- 
lytic in the absence of the SOI32&M Ag 

our calculation 
shows, the answer to this question is negative. 

One more comment on the spin and charge susceptibil- 



ities is in order: despite the fact that we considered only 
the Rashba SOI, all our results are applicable to sys- 
tems where the Dresselhaus SOI with coupling strength 
j3 takes place of Rashba SOI, i.e., (3 ^ 0, a — 0; in this 
case, the Rashba SOI should be simply replaced by the 
Dresselhaus SOI (a — > (3). 

Finally, we analyzed the nonanalytic dependence of 
the free energy, J 7 , in the presence of the SOI and at 
zero temperature beyond the Random Phase Approxi- 
mation (RPA). This analysis is important in the con- 
text of interacting helical Fermi liquids that have recently 
attracted considerable attention.— ~— In contrast to the 
RPA result^ 2 -, which predicts that the free energy scales 
with a as a 4 ln|a|, our result shows that the renormal- 
ization is stronger, namely, T oc C/ 2 |a| 3 C(J71n \a\), where 
C(x -> 1) ~ x 2 and C(x -> oo) ~ 1/x 2 . 
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Appendix A: Derivation of common integrals 



In this Appendix, we derive explicit expressions for some integrals of the Green's function which occur throughout 
the paper. 



1. "Quaternions" (I; mnr and J; mnr ) and a "triad" (Ii mn ) 



The first integral is a "quaternion" -a convolution of four Green's functions defined by Eq. (|2.7e[) . This convolution 
occurs in diagram 1, where it needs to be evaluated at small external and transferred momenta: q, q -C kp. To linear 
order in q and a, £k+ q + sa|k + q| = ej. + vpq cos 6k q + akp + o(q 2 , aq) with 9k q = Z(k, q). The same approximation 
holds for q with 9 k q = ^(k, q). Switching to polar coordinates and replacing kdk by m*dek, we reduce the integral to 



d6 kq f du)k f de k I ! 



2ir J 2ir J 2n iu>k ~ £fe — VpqcosOkq — Ictkp iuj^a — — makp 

x I \ (Al) 

i(uik + Q) — tk — VFqcosOkq — nakp iu>k — eu — rakp 



Integrating first over and then over w/., we obtain 



(2n) 2 J iQ — VFqcosOkq + (m — n)akp zf2 — VpqcosOkq + (r — n)akp 

(A2) 

id — vpqcos9kq + VFqcos9k q + (I — n)akp 
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Finally, the integral over 9k q gives 



2n (r — m)afci?[(Z — m)akp + Vpq cos 9k^][(l — r)akp + vpqcos9kq] 
(I — r)akp + vpq cos 9kq (I — m)akp + vpqcos9kq (r — m)akp 



\J vpq 2 + (17 + i(n - m)akp) 2 y/vpq 2 + (fi + i{n — r)akp) 2 \Jv 2 F q 2 + (17 — ivpqcos9kq + i(n — l)akp) 2 

(A3) 

Because to the overall term (r — m)akp in the denominator, the case r = m has to be treated specially. Taking the 
limit Jj mnm (0, 9 k q, q, q) = lim,.^ m I/„ mr (Sl, k q, q, q), one obtains 



J -lmnm (n,9 k q,q,q) 



mi 



27r [(Z — m)akp + vpqcos 9kq] 



1 



\/v' F q 2 + (fi — ivFqcos9kq + i{n — l)akp) 2 



v^q 2 + [17 + i(n — m)akp ][17 + ivpq cos 9kq + i(Z + u — 2m)aki 
[v 2 F q 2 + (fl + i{n - m)ak F ) 2 f /2 



Similarly, we obtain for another quaternion Ji mn r, defined by Eq. (|2. 15c() 
m*|17| 1 
imnr( , feq, 9, 3 2^ — ivpq cos 9kq -\- i(n — m)akp] [17 — ivpq cos 0fcg + i(r — ^)afci?] 

1 + 1 

\/v 2 F q 2 + (17 — ivpqcos9k q + i(r — m)akp) 2 \/vpq 2 + (17 — ivpqcos9k q + i(n — l)akp) 2 

Finally, we obtain for a convolution of three Green's functions-a " triad" -defined by Eq. (|2.15d[) 

m*\n\ 1 



Ilmn(£l,Okq,q,q) 



2ir vpq cos 9kq + (I — m)akp 
1 



1 



\Jv 2 F q 2 + (17 + i(n — tjakp — VFqcos9kq) 2 \/v 2 F q 2 + (17 + i(n — m)akp)' 2 



(A4) 



(A5) 



(A6) 



2. Integrals over bosonic variables 

There is a number of integrals over the bosonic frequency 17 and momentum q one encounters while calculating 
the spin susceptibility. The following strategy provides a convenient way of calculating all of them: (i) integrate over 
VFq for x £ [0,oo [, (ii) integrate over 17 by introducing a cut-off A-the low-energy physics proves to be independent 
of the choice of the cut-off, (iii) perform angular integration, which is trivial for the out-of-plane spin susceptibility 
and, in that case, can be performed at the very beginning. 

Again, it is convenient to treat the out-of-plane and in-plane components separately. 



a. Out-of-plane components 

As it was explained in the main text, the q dependence of \ zz f° r ? ^ 9q can be calculated perturbatively, by 
expanding in q/q a , where q a = 2m*|a|. In this section, we calculate only the leading term of this expansion obtained 
by setting q = 0. Corrections enter only quadratically in q and a more detailed calculation is necessary in order to 
show that. Later, in Appendix [B] we find the entire dependence of \ zz on <Z exactly, and show that this dependence 
is absent for q < q Q , which means that all terms of the expansion in q/q a vanish. For now, we focus on the q = case 
and evaluate the integrals in Eqs. (|2.8a[) and (|2.8b|) for xT 

^ fd 9 M r<Kk {I+ ___ +I +++)Il 



2?r / 2?r / 2tt 



o — 

ail I xdx — — — , + c.c. 



%-K 2 v F uk F ) J Vz 2 + 17 2 \_ ^Jx 2 + (17 + 2iak F ) 2 Vx 2 + 17 2 

/ m \ 2 f A inn2 , 17 2 ( m \ 2 \a\kF .... 

' ' ' dtm 2 h, 9 ■ — 5 = — — U— + ... (A7) 



\8Tr 2 VFCikF J j_a fl 2 + a 2 k F \4itvf J 67r 
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dn r do k g 

2tt / 2tt 



f^(/_ + _ + n_ + + / + _ + _n + _) 



m 



( 1 



dSX \ xdx 



y/x 2 + (n + 2iak F ) 2 

2iakp(n — 2iakp) 



Vx 2 + n 2 y?x 2 + (fi - 2iak F ) 2 [x 2 + 2iak F ) 2 fl 2 



l-K 2 VFOikp 



dflfl 2 In 



n 2 



n 2 



Y Hi 



i 2 kp \Attvf J 6n 



+ c.c. 



(A8) 



where . . . stands for non-universal, A-dependent terms and c.c. denotes the complex conjugate of the preceding 
expression. Substituting these results back into Eqs. (|2.8al) and (|2.8b[) . we obtain Eqs. (|2.9a[) and (I2.9b|) . Similarly, 
we obtain for the combination of triads in Eq. (|2.18[) for x¥' 



2 =- 



dn r do 



2tt / 2tt 



qdq 

Z7T 



-+++ I+-+I-+-) 
2 



m 



4:TT 2 VFakF 

/ \ 2 

/ m 

\4:Tr 2 VFOikF 



dft / xdx 



dflfl 2 In 



n 2 



yjx 2 + (n + 2iak F ) 2 Vx 2 + n 2 

2 ( m V 



n 2 + a 2 k 2 F 3tt \4ttvf) 



(A9) 



b. In-plane components 



We start with X\ x given by Eqs. (|2.11al) and (12. 1 lb[) . First, we notice that the quaternion structure of the first lines 
in Eqs. (|2.11al) and (|2.11bl) is the same as in the first lines of Eqs. (12.8a[) and ()2.8bl) for the out-of-plane component; 
the only difference is in the factor of sin 2 Okq- Since these expressions contain a, they can be evaluated at q = in the 
same way as the corresponding expressions in \ z \ were evaluated. At q — 0, the factor of sin 2 6% q just gives 1/2 of the 
corresponding contribution to \\ z ■ Next, we calculate explicitly the integrals in the second line of Eq. (|2.11a[) and in 
the third line of Eq. (|2.11bl) . These contributions contain an overall factor of q~ 2 and, therefore, one has to calculate 
the full dependence on q without expanding in q/q a . The part of the integrands that are odd in the angle drop out 
and, since d9f(i cos 9) — d6[f(i cos 0) + /(— icos#)]/2, all the formulas can be written in an explicitly real 
form. For the first of these two integrals we obtain (for brevity, we relabel O^q — > 9) 



f^fMfqJq 
J 2vr J 2vr J 2tt V 4 



A.ix 2 v 2 F q 
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— / dVL 
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n 2 
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ATT 2 Vpq 
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^x 2 + n 2 



c.c. = — 



A-K 2 v 2 F q 
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2 ,-A 



/ dn 



-A 



2 7T 



2tt 
dO 



xdx- 
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\Jx 2 + ft 2 \ y^X 2 + (ft + iVF<l cos ( 



dn I —n 2 in ( i 

2tt 
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cos 0' 



^Air 2 v 2 F q 
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v F q 2 
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/ tti \ VFq 
\4itvf J 97T 2 



(A10) 
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where, as before, . . . stands for non-universal, A-dependent terms. Notice that the SOI dropped out and, therefore, 
the equation above is valid for any ratio q/q a - The second integral reads as 



dfl f dOkq f qdq 
2tt J 2tt J ~2tt 

2 



cos 2 e(i ++ -+ii- + + + _n+_) 



m 



Air 2 v 2 F q 



— / dQ I .rd.r 
1 



n 2 



y/x 2 + (Q + 2iak F ) 2 



\J x 2 + (ft — 2iakp — ivpqcosd) 2 \/ x 2 + (57 — 2iakp + ivpq 



i£l 2 vpqcos9 



\x 2 + (fi + 2iak F ) 2 \< 



57 + 2iakp 



x 2 + + 2iak F y 



cos(9) 2 yjx 2 + (fi - 2iak F ) 2 
dQ 

2^'"V ' 4S7 2 



4<?/ J -A 



di m 2 I ^ln f 1 + ^J- COS 2 i 



f m \ v F q 



\<Ittv f ) 



9tt 2 



(All) 



which is the same result as in Eq. (|A10[) . The second line in Eq. (|2.11b[) gives the same result as the third one. 
Collecting all the results above, we arrive at Eqs. (|2.12a|) and (|2.12bj) . 

Finally, for the SOI-independent part of diagram 3, we find 



dfl f d6f.g f qdq 



kq 

2tt I ~2tT 



2tt 



cos 2 0[/ +++ (fi, 9 k g,q, g)/___(fi, 9 k g, qi -q) + /++_ (fi, 9 m , q, $)/__ +(fi, 9 k g, q, -q) + (q -> -q)\ 



■} F q 
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\J x 2 + (57 + iv F qt 



^x 2 + n 2 



i 



= -2 



^4ir 2 v F q 

( m \ 2 vpq 
~ \ 2nv F J 9TT 2 



d9 



dfifi 2 / — In 1 



2 ~2 9 

v F q cos 

fi 2- 



y/x 2 + (57 + 2iakp + iv F q cos 9) 2 y/x 2 + (57 + 2iak F ) 2 

2 „A 



= -2 



2n 2 v 2 



v F q 



dnn 2 in - ( l 



i 



4S7 2 



(A12) 



Appendix B: Full q dependence of the spin susceptibility 

In the main text and preceding Appendices we found \ zz a t zero external momentum. Here, we show how the full 
dependence of \ zz can be found using the q = part of diagram 1 in Fig. [2] as an example. 

We consider Eq. (|2.8a[) at finite q. The integral over bosonic variables reads as 



dtt f d9 
2tt J 2tt 

( m 
\4:Tr 2 v F 



qdq 



[I- 



+- 



2?r 

" dO f°° f°° 

27T./-OC Jo Vx 2 + 57 2 



2tt 
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+++ }u (n, q ) 
S7 2 
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i 



(2ak F + vpq cos 9) 2 
iS7(2a + vpqcos9) 



\Jx 2 + (57 — 2iakp — ivpqcos9) 2 \Jx 2 + fl 2 



( x 2 + Q2)3/2 



+ (a ->• -a) 



(Bl) 



where (a — > —a) stands for a preceding term with a reversed sign of a and, as before, we relabeled 9 k q — > 9. The last 
term in the parenthesis vanishes upon integration over either the angle (in the principal value sense) or the frequency 
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(it is odd in O), whereas the remainder yields 
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cdx 



n 2 
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\Jx 2 + Q 2 {2ak F + v F qcosl 
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I 6 [ 
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for q <q a , 
for q^> q a . 



(B2) 



We see that while xf z is independent of q for q < q a , for g Q 3> g Q it approaches the linear- in-g form found in Ref. [20| 
in the absence of the SOI. 

Another integral of this type occurs in the in-plane component, e.g., in the first line of Eq. (I2.11ap . The only 
difference compared to the out-of-plane part is an extra sin 2 9 factor. The q and f2 integrals are calculated in the 
same way while the angular integral is replaced by 



dn f d9 f qdq . 2 



2tt 



2tt 



n 



24 



/ d0sin 2 0|2afcF +«Fgcos0| = -^Re 
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2+1 ^ 
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1+| % 



for q <q a , 
for q^> q a - 



(B3) 



Appendix C: Logarithmic renormalization 

In this Appendix, we analyze renormalization of the out-of-plane component of the spin susceptibility in the Cooper 
channel for q <C q a - As an example, we consider diagram 1 at large momentum transfer to third order in the electron- 
electron interaction, see Fig. [9j The calculation is carried out most conveniently in the chiral basis as shown below 



xT q =2 kF = ~2(7 3 J Q J K J^J^J p J L Tr i G (- K + Q)G{-P + Q)G(-L + Q)]Tr[G(K + Q)a x G(K)G(P)G(L)G(K)a x ] 

= —4U 3 J — J j 7j~ J 2~ ^ ^tm;ns(0k,9p)T ns - vu (9 p , 9i)T vu .t r {9l,9k)^rl^k)^l m {^k)hmnr^stL UVl 



(CI) 



where 

F „ ,J0u.9J\ = T7(n. .sJk. s,Un. .sdk..<?o\ = 

is the scattering amplitude in the Cooper channel (k = — p) , a s t(9k) = (k, s| a 31 |k, i) = —i(se iek — te~ l9k )/2, and 



r siS2;S4S3 (0 fc , P ) = C/(p, s 3 |k, Sl )(p, S4 |k, s 2 ) = ^ ( 1 + SlS3 e^-^ ) ( 1 + S2S4 e^-^ ) (C2) 
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2^ 



du>i / deig u (L)g v (-L + Q) 



A 2 f £ln-4p- = L(a) for u = « 



= 7T ln /,. — T7i I 77. 7X77. I m = < m , a ^ _ r ,x r. ( C3 ) 



4tt (^ffcos^ + Cu-^afe^ + n 2 U^^ivni^W for « 



- r 



is the particle-particle (Cooper) propagator, evaluated on the Fermi surface at fixed direction of the fermionic mo- 
mentum 1. An additional factor of 2 in Eq. (|C1I) is related to the possibility of extracting the logarithmic contribution 
from either the integral over P or that over L. Note that each scattering amplitude depends on the difference of two 
angles, i.e., 9 p — 9\ = 9 pq — 9i q , such that the angle 9\ q is shared between the vertices and the function L uv . Moreover, 
due to the correlation of momenta, we have 9 P = 9k + ft and 9k — n/2 — 9kq. 
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-K + Q, t 




K + Q,l 

FIG. 9: Diagram 1 to third order in electron-electron interaction at large momentum transfer; here, the in-plane component is 
shown. 



Upon summation over the Rashba indices, the integration over 9i q is readily carried out in all u = v terms, whereas 
the u = —v terms require more a careful treatment. Due to the dependence of the scattering amplitudes on 9i q , L uv 
enters multiplied by a either constant, or by sin26'/ g , or else by cos 29i q 
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d6 lq 
2ir 



1 

sin 29 iq 
cos 26l q 



Obviously, only the first choice leads to logarithmic renormalization. Keeping only this choice for u 



In 




(C4) 



we obtain 



XX 

X.X,q=2ki 



= - u 6 — / — 



sin 2 6» fe( j(/ + _ H 
sin 2 9k q (I+- 
3 sin 2 9 kq (I+- 



d6 



kq 



2tt 



qdq( 
2tt \ 



_n+_ 



3 sin 2 9 kq (I+ + 7- + ++)n 

n_+) + cos 2 fe g(i ++ _ + n_ + + 



3 cos 2 e kq (i +++ 
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i )n (C5) 
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-)n 



cos 



h q {I. 



++-H 



-)n 



/_+_+n_ + ) + 3 cos 2 fe9 -(/++_ + iL 



n. 



iln- 



A 



(C6) 



The first two lines in Eq. (|C6|) contain an Q and q- independent logarithmic factor. Integrations over q, 9 kq , and f2 
in these lines produce terms which scale either as q or as |a| , thus these two lines generate terms of the type gin \ a\ 
and \a\ In |a| . Next, we note that some combinations of quaternions and polarizations bubbles in these two lines, when 
integrated over q, 9 kq , and fl, produce a q term while others produce an |a| term. Namely, combinations (/ ++++ + 



-)n and I ++ h n_ + +7 | II + _ produce q , while (7 H +I_ 



+++ 



)n and I 



+-+ 



+-+ 



II |_ produce 



a\. To extract the leading logarithmic dependence, we split the Q and q-dependent logarithmic factor into two parts as 

^F*i Lin ^ tttVi^ivi i+ TYml + TT\licic +V>Q (inTYiKinofiATio t~\f +tio fir-of tuno and oo In kj^ \ct \ 



In 



; +ln — when it multiplies the combinations of the first type, and as In 

M+y'v 2 F q 2 +n 2 JV ' 



+ In , 4 , when 



it multiples the combinations of the second type. The f2- and q-dependent remainders do not produce main logarithms 
because the internal scales of ft and q are set either by q or by |a| for the first and second types, correspondingly. 

Therefore, the only main logarithms we have are either In or In fcF ^ a ■ Collecting all the contributions, we finally 
obtain 
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v F q 



(C7) 



where in the last line we retained only leading logarithms renormalizing each of the two terms in of the second-order 
result. Thus we see that each energy scale, i.e., vpq and VFq a , isrenormalized by itself. 
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FIG. 10: Left:A skeleton diagram for the free energy in the presence of the Cooper renormalization; f is a renormalized Cooper 
vertex. Middle: The effective scattering amplitude ri 4 s 2;S3S4 (k, k'; p, p') in the chiral basis. Right: A generic n-th order ladder 
diagram in the Cooper channel, ri™s 2;S3S4 (k, — k; p, — p). 



Appendix D: Nonanalytic dependence of the free energy as a function of spin-orbit coupling 

In a number of recent papers,— ~— the properties of interacting helical Fermi liquids were analyzed from a general 
point of view. In particular, Chesi and Giuliani^ have shown that an equilibrium value of helical imbalance 



SN = 



N + -N_ 
N + + N_' 



(Dl) 



where N± is the number of electrons in the ± Rashba subbands, is not affected to any order in the electron-electron 
interaction and to first order in Rashba SOI. Mathematically, this statement is equivalent to the notion that, for small 
SN and a, the ground state energy of the system T can be written as 8 = A(SN — 2ma/fcf) 2 , so that the minimum 
value of T corresponds to the non-interacting value of SN. The analysis of Ref. was based on the assumption that 
JF is an analytic function of a, at least to order a 2 . In a related paper, Chesi and Giuliani^ analyzed the dependence 
of J- on SN within the Random-Phase Approximation (RPA) for a Coulomb interaction and found a non-analytic 
^ 4 In \SN\ term. 

In this Appendix, we analyze the non-analytic dependence of J 7 on a by going beyond the RPA. [For small a, 
there is no need to consider the dependencies of J 7 on a and SN separately, as the shift in the equilibrium value of 
SN due to the electron-electron interaction can be found perturbatively.] To this end, we derive the free energy at 
q = T = 0-equal, therefore, to the ground state energy-following the method of Ref. which includes renormalization 
in the Cooper channel to all orders in the interaction. 

The free energy is given by the skeleton diagram in Fig. 1101 



*^ zz — 



,(k F ,-k F ;-k F ,k F )f s 



l(— k F , k F ; k F , — k F )n siS2 iT S3S4 , 



(D2) 



where a particle-hole bubble is given by Eq. (|2.7fj) and f SlS2 - S3S4 (k F , — k F ; — k F , k F ) is a scattering amplitude renor- 
malized in the Cooper channel. To first order in electron-electron interaction U, T SlS2 - S3Si is given by Eq. (|C2|) . 
It is convenient to decompose the renormalized amplitude into s, p, and d channels as 



r« 2;S3S4 (k,-k;p,-p)(L) = ^ lS2;S3S4 (L) + ^ lS2;S3S4 (L) e ^-^ +WW 3S4 (L) e 2 ^-H 



(D3) 



where the bare values of the corresponding harmonics are U SlS2 - S3Sll (0) = U2fc F /2, V SlS2 - S3Si (0) = W2/c F (siS3 + S2S4V2, 
and W S1S2;S3S4 (0) = U2k F SiS2S^,s i /2. The s,p,d harmonics of f were shown in Ref. to obey a system of decoupled 
Renormalization Group (RG) equations: 
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where the RG variable is defined as 
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(D4) 
(D5) 
(D6) 

(D7) 



and the initial conditions were specified above. Solving these equations, we obtain U SlS2 - S3Si (L) — u/[2(l + uL)], 
V SS ;± S±S (L) = ±u, V SlS2S3Sll (L) = u(sis 3 +s 2 S4)/(1+2mL) for the remaining Sj's, and W SlS2 - S3Si (L) = usis 2 s 3 s 4 /[2(l+ 
uL)], with u = U2k F - Combining the solution in the Cooper channel, we find 



f s ± S ;± ss (k F , -k F ; -k F , k F ) 



u 



1 + uL 

~ u 
T s _ s .± STS (k F , -k F ; -k F , k F ) = - i — =p 



1 + uL 1 + 2uL 

and zero for the remaining cases. 

Substituting the RG amplitudes into Eq. (|D2j) and summing over the Rashba indices, we arrive at 
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(D10) 



The terms proportional to I1q are divergent and scale with the upper cut-off A, thus they can be dropped as we are 
interested only in the low energy sector. Making use of the following integrals J qdq(H^ II |_ — IIq) = and 



J dfin 2 J dqq(Il 2 + _ + Ii 2 _ + - 2U 2 ) = J dtt%- In 



a 2 



r fl 2 + 4a 2 k 2 F 



16tt 
34 



a\ 3 k 3 F + 0(A), 



(Dll) 



we obtain the final result 
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Note that T is non-zero starting only from the fourth order in u: 
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\a\ 3 k 3 F j 2 fVA^i 
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(D13) 



Apart from the logarithmic factor, a cubic dependence of T on \a\ is in line with a general power-counting argument 24 
which states that the non-analytic dependence of the free energy in 2D is cubic in the relevant energy scale. A cubic 
dependence of J 7 on a implies that the shift in SN scales as a 2 C(L), where C(L) is a function describing logarithmic 
renormalization in Eq. (|D12[) . This is to be contrasted with an a 3 In a scaling predicted within the RPA42 
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